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Abstract — This work presents blind constrained constant mod- 
ulus (CCM) adaptive algorithms based on the set-membership 
filtering (SMF) concept and incorporates dynamic bounds for 
interference suppression applications. We develop stochastic 
gradient and recursive least squares type algorithms based on 
the CCM design criterion in accordance with the specifications 
of the SMF concept. We also propose a blind framework that 
includes channel and amplitude estimators that take into account 
parameter estimation dependency, multiple access interference 
(MAI) and inter-symbol interference (ISI) to address the impor- 
tant issue of bound specification in multiuser communications. A 
convergence and tracking analysis of the proposed algorithms is 
carried out along with the development of analytical expressions 
to predict their performance. Simulations for a number of 
scenarios of interest with a DS-CDMA system show that the 
proposed algorithms outperform previously reported techniques 
with a smaller number of parameter updates and a reduced risk 
of overbounding or underbounding. 

Index Terms — Interference suppression, blind adaptive estima- 
tion, set-membership estimation, spread spectrum systems. 



I. Introduction 

Set-membership filtering (SMF) E), Js], JS], iU is a class 
of recursive estimation algorithms that, on the basis of a 
pre-determined error bound, seeks a set of parameters that 
yield bounded filter output errors. These algorithms have 
been applied to a variety of applications including adaptive 
equalization lO and multi-access interference suppression |6l, 
EJ. The SMF algorithms are able to combat conflicting 
requirements such as fast convergence and low misadjustment 
by introducing a modification on the objective function. These 
algorithms exhibit reduced complexity due to data-selective 
updates, which involve two steps: a) information evaluation 
and b) update of parameter estimates. If the filter update 
does not occur frequently and the information evaluation 
does not involve much computational complexity, the overall 
complexity can be significantly reduced. 

Adaptive SMF algorithms usually achieve good convergence 
and tracking performance due to the use of an adaptive step 
size or an adaptive forgetting factor for each update. This 
translates into reduced complexity due to the data selective 
updating. However, the performance of SMF techniques de- 
pends on the error-bound specification, which is very difficult 
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to obtain in practice due to the lack of knowledge of the 
environment and its dynamics. In wireless networks character- 
ized by non-stationary environments, where users often enter 
and exit the system, it is very difficult to choose an error 
bound and the risk of overbounding (when the error bound is 
larger than the actual one) and underbounding (when the error 
bound is smaller than the actual one) is significantly increased, 
leading to a performance degradation. In addition, when the 
measured noise in the system is time-varying and the multiple 
access interference (MAI) and the intersymbol interference 
(ISI) encountered by a receiver in a wireless network are 
highly dynamic, the selection of an error-bound is further 
complicated. This is especially relevant for low-complexity 
estimation problems encountered in applications that include 
mobile units and wireless sensor networks Jg], llTOl . where the 
sensors have limited signal processing capabilities and power 
consumption is of central importance. These problems suggest 
the deployment of mechanisms to automatically adjust the 
error bound in order to guarantee good performance and a 
low update rate (UR). 

A. Prior and Related Work 

In this context, blind methods are appealing because they 
can alleviate the need for training sequences or pilots, thereby 
increasing the throughput and efficiency of wireless net- 
works. In particular, blind algorithms based on constrained 
optimization techniques are important in several areas of 
signal processing and communications such as beamforming 
and interference suppression. The constrained optimization 
required in these applications deals with linear constraints 
that correspond to prior knowledge of certain parameters 
such as the direction of arrival (DoA) of user signals in 
antenna array processing ||231 and the signature sequence of 
the desired signal in DS-CDMA systems |l24|, ESI. There- 
fore, linear signal models and constraints can be used to 
describe various wireless communications systems including 
multi-input multi-output (MIMO) and orthogonal frequency- 
division multiplexing (OFDM) systems. For instance, linear 
constraints that incorporate the knowledge of user signatures 
of a DS-CDMA system can also be used to exploit the 
knowledge of the spatial signatures of MIMO systems. A 
number of blind algorithms with different trade-offs between 
performance and complexity have been reported in the last 
decades ll24l - ll36l . Local scattering, synchronization and es- 
timation errors, imperfectly calibrated arrays and imprecisely 
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known wave field propagation conditions are typical sources 
of uncertainties that can lead to a performance degradation 
of blind algorithms. The literature indicates that the CCM- 
based algorithms Il30l - ll36l have a superior performance to 
algorithms based on the constrained minimum variance (CMV) 
criterion ll24l - ||271 . The CCM-based algorithms exploit prior 
knowledge about the constant modulus property of signals like 
M-PSK, which results in an improved performance over CMV- 
type techniques and a performance that is very close to the 
linear minimum mean-square error (MMSE) training-based 
techniques. Moreover, the CCM-type algorithms are robust 
against errors in the effective signature sequence required for 
blind parameter estimation, which prevents a severe perfor- 
mance degradation in the presence of uncertainties. These 
features make CCM-type algorithms excellent candidates for 
interference suppression in wireless networks. The need for 
the adjustment of parameters (step size, forgetting factor) and 
the computational complexity of CCM-based techniques calls 
for approaches like SMF. Prior work on SMF blind algorithms 
for interference suppression is very limited ||7], ||8]. The use 
of time-varying bounds is also restricted to applications where 
one assumes that the "true" error bound is constant llT6l and 
to the parameter-dependent error bound recently proposed in 
El, ifTsl . The time-varying bound techniques so far reported 
are not bUnd and do not exploit these mechanisms for channel 
and parameter estimation. 



B. Contributions 

In this work, we propose set-membership blind adaptive 
constrained algorithms based on the CCM criterion. Prelim- 
inary results have been reported in ll22l . In particular, we 
derive stochastic gradient (SG) and recursive least squares 
(RLS)-type CCM algorithms designed in accordance with the 
specifications of the SMF concept. The second contribution is 
a low-complexity blind framework for parameter estimation, 
and tracking of parameter evolution, MAI and ISI that relies 
on simple estimation techniques and employs the proposed 
set-membership CCM algorithms with a time-varying bound. 
The third contribution of this work is an analysis of the 
optimization problem that gives rise to the proposed algorithms 
along with a mean-squared error (MSB) convergence and 
tracking analysis using the energy conservation approach ll39l 
for predicting the performance of the algorithms. A simulation 
study considers an interference suppression application to DS- 
CDMA systems, which compares the performance of the 
proposed and existing algorithms, and discusses the main 
features of the algorithms. 

This paper is structured as follows. Section II describes 
the DS-CDMA system model and briefly reviews linearly 
constrained receivers. Section III introduces the SM blind al- 
gorithms with time-varying bound. Section IV proposes blind 
parameter dependent and interference dependent bounds along 
with blind channel and amplitude estimation algorithms. Sec- 
tion V is dedicated to the analysis of the proposed algorithms. 
Section VI is devoted to the presentation and discussion of 
numerical results, while Section VII gives the conclusions. 



ii. ds-cdma system model and linearly 
Constrained Receivers 

Let us consider the uplink of a symbol synchronous DS- 
CDMA system with K users, N chips per symbol and Lp 
propagation paths. A synchronous model is assumed for sim- 
plicity since it captures most of the features of more realistic 
asynchronous models with small to moderate delay spreads. 
The modulation is assumed to have constant modulus. Let 
us assume that the signal has been demodulated at the base 
station, the channel is constant during each symbol and the 
receiver is perfectly synchronized with the main channel path. 
The received signal after filtering by a chip-pulse matched 
filter and sampled at chip rate yields an il/-dimensional 
received vector at time i 
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where M = N + Lp — 1, n[i] ~ [ni[i] ... 
is the complex Gaussian noise vector with zero mean and 
i?[n[i]n^[i]] — a^I independent and identically distributed 
samples, {.)'^ and (.)^ denote transpose and Hermitian trans- 
pose, respectively, and E[] stands for expected value. The user 
symbols are denoted by hk[i], the amplitude of user k is Afc[i], 
and is the intersymbol interference (ISI) for user k. The 
signature of user k is represented by Ck = [ca;(1) . • • Ck{N)]^ , 
the MxLp constraint matrix Ck that contains one-chip shifted 
versions of the signature sequence for user k and the Lp x 1 
vector hk [i] with the multipath components are described by 



Ck = 
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, hk \i] = 



hk,Lp-i\i] 



(2) 

The MAI comes from the non-orthogonality between the 
received signature sequences. The ISI originates from the 
multipath propagation effects of the channel, depends on the 
length of the channel response and how it is related to the 
length of the chip sequence. We define Lg as the ISI span, 
i.e., the number of symbols affected by the channel. For 
Lp = 1, Ls = 1 (no ISI), for I < Lp < N,Ls = 3, 
for N < Lp < 2N, Lg = 5, and so on. At time instant i 
we will have ISI coming not only from the previous time 
instants but also from the next symbols. The linear model 
in ^ can be used to represent other wireless communica- 
tions systems including multi-input multi-output (MIMO) and 
orthogonal frequency-division multiplexing (OFDM) systems. 
For example, the user signatures of a DS-CDMA system are 
equivalent to the spatial signatures of a MIMO system. 

In order to describe the design of linearly constrained 
receivers, we consider the received vector r[i], the AI x Lp 
constraint matrix Ck that contains one-chip shifted versions of 
the signature sequence for user k and the Lp x 1 vector hk [i] = 
[hkfii'i] ■■■ hk^L -i[t]V with the multipath components to 
be estimated. The CCM linear receiver design is equivalent 
to determining an FIR filter Wk[{\ with M coefficients that 
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provide an estimate of the desired symbol as follows 

Zk[i]=wj^[i]r[i], (3) 

where the detected symbol is given by bk[i] = Q{wj^[i]r[i]), 
where Q( ) is a function that performs the detection according 
to the constellation employed. 

The design of the receive filter Wk[i] is based on the 
optimization of the CM cost function 



JcMiw,[i])=E\{\wj^[i]r[ir-lf 



(4) 



subject to the constraints given by = where 

Pk[i] = Ckhk[i] is the effective signature vector that 
corresponds to the convolution between the original signature 
sequence and the channel gains, and is a constant to 
ensure the convexity of the optimization problem as will be 
discussed later on. This approach assumes the knowledge 
of the channel. However, when multipath is present these 
parameters are unknown and time-varying, requiring channel 
estimation. The CCM receive filter expression that iteratively 
solves the constrained optimization problem in (|4]i is given by 



Wk[i + 1] =-Rfc^M dk\i]- 



1 2 



(5) 



where Zk[i\ = w^[i\r[i], Rk[i] = E[\zk\i]\^r[i\r" 
dk[i] = E[zl[i]r[i]]. A detailed derivation of the CCM 
estimation approach can be found in ll33l . ll34l . 1361 . It should 
be remarked that the expression in (|5) is a function of previous 
values of the filter Wk [i] and therefore must be iterated in order 
to reach a solution. In addition to this, the iterative method in 
([5]l assumes the knowledge of the channel parameters. Since 
there is a large number of applications that have to deal with 
unknown multipath propagation, it is also important to be able 
to blindly estimate the multipath components. 

In order to blindly estimate the channel, a designer can 
adopt the blind channel estimation procedure based on the sub- 
space approach reported in ||261 . 11411 and which is described 
by 

hk[^ - argmin h^[i]C^ R-^[i]Ckhk[i] (6) 

subject to \\hk[i]\\ = 1, where R[i] = £'[r[z]r^[i]]. The 
solution is the eigenvector of the L„ x Lp matrix corresponding 
to the minimum eigenvalue of Cj. R^^[i]Ck obtained by an 
eigenvalue decomposition (EVD). Here, we use Rk[i] in lieu 
of R[i] to avoid the estimation of both R[i] and Rk[i], which 
shows no performance loss as reported and investigated in 

1331, m, M- 

III. Set-Membership Blind Adaptive Constrained 
Algorithms with Time Varying Error Bounds 

In this section, we describe an adaptive filtering framework 
that combines the set-membership (SM) concept with blind 
constrained algorithms based on the CCM design. We also 
introduce simple time-varying error bounds to take into consid- 
eration the evolution of the receive filter and the MAI and ISI 



effects. In SM filtering ||2l, the parameter vector Wk [i] for user 
k is designed to achieve a specified bound on the magnitude of 
an estimated quantity Zk [i] ■ As a result of this constraint, the 
SM bUnd adaptive algorithm will only perform filter updates 
for certain data. Let 8 [i] represent the set containing all Wk [i] 
that yield an estimation quantity upper bounded in magnitude 
by a time-varying error bound jkii]- Thus, we can write 

0W= n WeC^^:|zfe[z]|<7fe[z]}, (7) 

{r[i])eS 

where r[i] is the observation vector, S is the set of all 
possible data pairs (6/c[i], r[i]) and the set 8[j] is referred 
to as the feasibiUty set, and any point in it is a valid 
estimate Zk[i] = wj^[i]r[i]. Since it is not practical to predict 
all data pairs, adaptive methods work with the membership 
set V'j = 0™=! ^™ provided by the observations, where 
Tim = {wk G C^' : \zk[i]\ < 7fc[*]}- In order to devise an 
effective SM algorithm, the bound jk [i] must be appropriately 
chosen. Due to the time-varying nature of many practical 
environments, this bound should also be adaptive and capable 
of estimating certain characteristics for the SM estimation 
technique. We devise SM-CCM algorithms equipped with 
variable step sizes and forgetting factors that are able to au- 
tomatically tune to different situations, which is an advantage 
over methods operating with fixed parameters in time-varying 
scenarios. In what follows, we derive SM-CCM algorithms 
for blind parameter estimation that assume time-varying error 
bounds. 

A. Set-Membership CCM Stochastic Gradient-Type Algorithm 

Here we develop the set-membership CCM stochastic 
gradient-type (SM-CCM-SG) algorithm. The basic idea is to 
devise a gradient descent strategy to compute a parameter 
vector Wk for user k that minimizes the instantaneous CM 
cost function subject to certain constraints, which requires for 
adaptation that the square of the error [i] exceeds a specified 
error bound 7|[i]. Mathematically, the proposed SM-CCM-SG 
algorithm solves the following optimization problem 

minimize JcM{wk[i]) (|if;f [i]r-[i]|^ - 1)^ = ei[i] 
subject to io^[i]pj,[i] ~ V (8) 
whenever e^[z] > 7^[i] 

This problem can be solved using the method of Lagrange 
multipliers using the equality constraint e^.[*] = 7fcW llll: 

L{wk^\,Kk) = (|u;f 
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i^kKPk [i'\Wk[i\ 



(9) 

where Kk is a Lagrange multiplier Computing the gradient 
terms of (|9]l and equating them to zero, we obtain 

'^■u,i[i\C{wk[i],Kk) = 2{ek[i]r[i]zl[i]) + Pk[i]Kk = 



\J,,'C{wk[i],Kk) = Wf, [i]Pk[i\ 



V 



0. 



We employ a gradient descent rule to solve for the above 
equations. Using the first equation of ( fTOl i. we have 

Wk\i^\\ =Wk\^ - pk{e.kW-l\i\r\i\ +Pfe[i]Kfe) (11) 
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Using the second equation of ( fTOl i. we obtain the exponentially-weighted CM cost function subject to con- 

„ ^traints that require that the squared error of the filter exceed 

Wk[i+1] = nkmwk[^-fi,e,[t]zl[i]r[i])+,.p,[i]{p, [i]p,[i]) ^ ^pg^jgg^ time-varying error bound 7^^. Mathematically, 

(12) 

2 



where is the effective step size, ek[i] = \zk[i 
1 is the error signal for user k, and Tlkii] ^ I 
Pj^[i]{pj^[i]Pi^[i])^^pjf[i] is a projection matrix that ensures 
the constraint and / is an identity matrix. By imposing the 
condition to update whenever e^[i] > 7^[z] we arrive at the 
set of all Wk that satisfy 



the proposed SM-CCM-RLS algorithm solves the optimization 
problem 



(13) 



minimize jQfj{wk[i]) = A 

subject to — v 

whenever e\{i\ > 



i-lr 



(21) 



It can be verified that the set above is non-convex and 
comprises two parallel hyper-strips in the parameter space. 
From the above conditions we consider two cases: i) \wj^[i + 
l]r[i]\ < v/l + 7fcM andii) \wj^[i + l]r[i]\ > ^1 - 7fcW- By 
substituting the recursion obtained in (fT2t into |t(;^[i + l]r[i]|, 
we have 



where \k[i] is a time-varying forgetting factor Similarly to 
the case of the SM-CCM-SG algorithm, this problem can be 
solved with the method of Lagrange multipliers [JJ along with 
the use of the condition e^[i] > 7^[i] to save computations 



\w"\i- 



L'^1]t-WI = \zk[i]~^^k[i]ek[i]zk[^]r"[^]Uk[^]r[i]\. (14) 
Using the above expression we have for case i): 

\zk\i] - Hkmek[i\zk[i]r" [i\Uk\i]r[{\\ = v^l + 7fc [{ 
which leads to 



Ciwk[^,ek) = J2X[7'm\wfmi]\'-lf 



1=1 



(22) 



[elip^l^Wkli]^,^)], 



(15) 



where is a Lagrange multiplier Computing the gradient 
terms of the Lagrangian in (l22l i and equating them to zero, 
we obtain 



1 



Zk\l 



ek\i]r^\i]Ilk\i]r[i\ 



(16) 



Using ( fT4l i we have for case ii): 

\zk\i] - iJ.k[i\ek[i\zk[i]r" \i]Uk\i]r[{\\ = \J\ 
which results in the following 

- Vi-7fcH\ 1 



V„.[,]/:(ti;,H,6,.) = Y.\r'mzkW - \){rW^{^Wk{i\+PkWk 
1=1 

Ve-£(wfc[i],efc) = HpfcH - 1/ = 0. 
(17) Solving for the above equations, we have 



(23) 



Wk\i - 



dkV]~(p^mk 



^^k[^\ = 



\Zk\l\ 



ek[i]r"[i]U.k[i]r[i] 



(18) 



v]PkW 



(24) 



The resulting SM-CCM-SG algorithm is described by 



where Zk[i] = w\}[i]r[i], Rk[i] = 

T:i=i>^k'mzk[i]?r[i]r"[ii dk[i] = T:i=i>^k'mmi]- 

Wk[i+l] ^Ilk[i]{wk[i]~lJLk[i]ek[i]zl[i]r[i])+vpj^[i]{^ we need to compute dk\i\ and this is performed by 

(19) the following recursion 



where 



dk\i\^dk[i-l\ + \k\i\z*k\i\r\i\. 



(25) 



\zk[i\\ 



e-k [i]r"[i]u.k[i]r[i] > \/(l + At ftis poiut, wc need to compute R^. [i] efficiently and this 

is done by applying the matrix inversion lemma [IJ, which 



■]rH[i\nk[i\r[. 



if\zkm<Vi^ 

otherwise 

(20) 

The SM-CCM-SG algorithm described in (O-dlOli requires 

M + UR(57\/ - 1) additions and M + 1 + UR(5A/ + 4) 



R, 



R, 



l + AfeH|zfcH|V^H^fc'[i-l]rK 



(26) 

The last step of the SM-CCM-RLS algorithm is the use of the 
multiplications per received symbol, where UR is the update condition e|[i] > 7^[i] to save computations and to adjust the 
rate. optimal \k[i\- In order to adjust \k[i], the authors in 16] have 

advocated a strategy that yields bounding ellipsoids that lead 
to a simple innovation check with linear complexity, which 
considers the cost function 



B. Set-Membership CCM RLS-Type Algorithm 

In this part, we derive the set-membership CCM recursive 
least-squares-type (SM-CCM-RLS) algorithm. The idea is 
to devise a least-squares method to calculate a parameter 
vector Wk used at the receiver for user k that minimizes 



C 



AfcW 



Xk\i] 



em 



l + Xkli]r"[i]R^ [i-l]rli] 



1 

(27) 
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The maximization of the cost function in dZTT i leads to the 
innovation check of the proposed SM-CCM-RLS algorithm: 



r«[i]H^ [i-l]r[: 






1 



otherwise. 



(28) 

The SM-CCM-RLS algorithm described in (|24li-(l26ll and ^ 
requires M + IJR{4AP + 3M) additions and M + 2 + 
UR(4M^ + 6M + 3) multiplications per received symbol. It is 
worth mentioning that the computational savings can be quite 
substantial if the algorithm operates with a low UR. 

IV. Blind Parameter Estimation and Time-varying 
Bounds 

This section presents a blind framework employed to com- 
pute time-varying error bounds •jk [i] based on parameter and 
interference dependency. The proposed blind framework is 
an extension of the approach in 1(2 1 il that computes time- 
varying error bounds, and performs interference estimation 
and tracking. In contrast to ET\ that considers training -based 
recursions, a blind procedure for estimating MAI and ISI 
power levels is presented with a set-membership blind channel 
estimator and a blind amplitude estimator, which are employed 
in the adaptive error bound for the SM adaptive algorithms. 

A. Parameter Dependent Bound 

Here, we describe a parameter dependent bound (PDB), 
that is similar to the one proposed in [|T31 and considers 
the evolution of the parameter vector Wk[i] for the desired 
user (user k). The PDB recursion computes a bound for SM 
adaptive algorithms and is described by 



7fc[i + !] = (!- ;3)7fcH + I3yfa\\wk\i\ 



(29) 



where /3 is a forgetting factor that should be adjusted to 
ensure an appropriate time-averaged estimate of the evolution 
of the power of the parameter vector iffc[i]. The quantity 
Q;||iOfc[i]|p(T^[i] is the variance of the inner product of Wk[i] 
with n[i] which provides information on the evolution of 
K;fe[i], where a is a tuning parameter and cr'^[i] is an estimate 
of the noise power This kind of recursion helps avoiding too 
high or low values of the squared norm of Wk \i] and provides a 
smoother evolution of its trajectory for use in the time-varying 
bound. The noise power at the receiver should be estimated 
via a time average recursion. In this work, we will assume 
that it is known at the receiver 

B. Parameter and Interference Dependent Bound 

In this part, we develop a blind interference estimation and 
tracking procedure to be combined with a parameter dependent 
bound and incorporated into a time-varying error bound for 
SM recursions. The MAI and ISI power estimation scheme, 
outlined in Fig. [1] employs both the RAKE receiver and the 
linear receiver described in (3) for subtracting the desired user 
signal from r[i] and estimating MAI and ISI power levels. 
With the aid of adaptive algorithms, we design the linear 



receiver, estimate the channel modeled as an FIR filter for the 
RAKE receiver and obtain the detected symbol bk[i], which is 
combined with an amplitude estimate [i] for subtracting the 
desired signal from the output Xk[i\ of the RAKE. Then, the 
difference dk[{\ between the desired signal and Xk[i] is used 
to estimate the MAI and ISI power 



Blind linear 
receiver Wfc[i] 



Zk[{\ 



Q{-) 



Blind channel 
estimation 



hfc[i] 



RAKE receiver 

ffcW 




Blind amplitude 
estimation 



MAI and ISI 




power estimation 


► 



Fig. 1. Block diagram of tlie proposed blind scheme. 



C. Blind Interference Estimation and Tracking 

Let us consider the RAKE receiver with perfect channel 
knowledge, whose parameter vector /^[j] = Ckhk[i\ for user 
k (desired one) corresponds to the effective signature sequence 
at the receiver, i.e. Ck[i\ = Ckhk[i]- The output of the RAKE 
receiver is given by 



Xk [I 



fk b 



Ak w bk w /f m [*] + E Aj [^b, [z] /f H c, [i 



desired signal 



K 

E- 



MAI 



ISI 



(30) 



where f^[i]ck[i] = Pk[i] and f"[i\cj[i\ = pi,j[i] for j ^ 
1. The symbol pk represents the cross-correlation (or inner 
product) between the effective signature and the RAKE with 
perfect channel estimates. The symbol represents the 

cross-correlation between the RAKE receiver and the effective 
signature of user j. The second-order statistics of the output 
of the RAKE in ( f30l ) are described by 

E[wm= Ai[i]pi[^E[\b,m 



K 



K 



Y,A][mbAmm. 



H~ ~H 



1=1 
3^k 



(31) 
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From the previous development, we can identify the sum 
of the power levels of MAI, ISI and noise terms from the 
second-order statistics. Our approach is to obtain instantaneous 
estimates of the MAI, the ISI and the noise from the output 
of a RAKE receiver, subtract the detected symbol in Q from 
this output (using the more reliable linear multiuser receiver 
{wk[i])) and to track the interference (MAI + ISI + noise) 
power as shown in Fig. 1. Let us define the difference between 
the output of the RAKE receiver and the detected symbol for 
user 1: 

K 



Xk[i\ = Ak[i]bk[i\Pk[i], where Pk[i] = Ckhk[i]. Let us now 
perform an eigen-decomposition on R 



K 



k=l 



= [(t>s 0n] 



A,, + cr2/ 
cr2j 



(36) 



dk [i] 



xk m 



where 0^ and ^„ are the signal and noise subspaces, respec- 
tively. Since 0^ and 0„ are orthogonal, we have the condition 

cl)^Xk[^ = 4,^Ak[i]bk[i]Pk[i] = ct,^Ak[i]bk[^Ckhk[i] = 0. 
Ak[j]bk[i] « 5] HSfcW +/f Mr7fc[t|ence, we have 

r = h^[i\ c^bi[i]Aim,,<t>nM^bk[t]Ck hk[i 



k=2 



ISI 



MAI 



(32) 



By taking expectations on and taking into account the 

assumption that MAI, ISI and noise are uncorrected we have: 

K 



k=2 



(33) 

where the above equation represents the interference power. 
Based on time averages of the instantaneous values of the 
interference power, we consider the following algorithm to 
estimate and track i?[|fife[i]p] 



v[t + l] = {1 - I3)v[i] + I3\dk[i] 



(34) 



where /3 is a forgetting factor To incorporate parameter 
dependency and interference power for computing a more 
effective bound as an alternative to replace (|29t . we employ 
the parameter and interference dependent bound (PIDB) 

7fc[z + l] = (1 - /?)7feM v^[i] + ^a\\wk\\^a^[i]) , 

(35) 

where v[i] is the estimated interference power in the 
multiuser system and r is a weighting parameter that must 
be set. The equations in (|34| | and dSSl l are time-averaged 
recursions that are aimed at tracking the quantities |(iA;[i]P 
and {\/Tv'^[i] + a\\w kW^ cr'^ [«]), respectively. The equations 
in ( |34] | and ( |35] ) also avoid undesirable too high or low 
instantaneous values which may lead to inappropriate time- 
vaiying bound 7fc [i] . 



D. Blind Channel and Amplitude Estimation 

Let us now present a set-membership blind channel esti- 
mation (SM-BCE) algorithm to design the RAKE and linear 
receivers. Consider the constraint matrix Ck that contains one- 
chip shifted versions of the signature sequence for user k 
defined in (2) and the assumption that the symbols bk[i] are 
independent and identically distributed (i.i.d), and statistically 
independent from the symbols of the other users. Consider the 
covariance matrix R = [r[i]r^[i\\ and the transmitted signal 



(37) 



The above relation allows to blindly estimate the channel 
hk\i]- To this end, we need to compute the eigenvector 
corresponding to the smallest eigenvalue of Tfe. It turns out 



that we can use the fact that \Ymp^oo{R/ c"^) ^ 



't>n4>n 

|36l and, in practice, it suffices to use p = 1 or 2. Therefore, 
to blindly estimate the channel of user k in the DS-CDMA 
system we need to solve the optimization problem 

= argmin/if [z]Tfc/ifc[i], subject to ||ft.fc[*]ll = 1^ 

(38) 

In order to solve (|38] | efficiently, we rely on the SM estimation 
strategy and a variant of the power method 13 3 i that uses a 
simple shift is adopted to yield the SM-BCE 



hk[i] = (/ - Tk[i]T k[i])hk[i ~ 1], 



(39) 



where Tk[i] = l/ir[Tfc[i]] and hk[i} ^ hk[i]l\\hk[i\\\ to 
normalize the channel. The quantity Tfe[i] is estimated by 

tfe[i] = tfe[z - 1] + \k\i]C^Al\i]Pl[i]Ak[i]Ck, (40) 

where \k[i] is a variable forgetting factor that is obtained 
by ( |28] | and -Pfe[i] is computed according to ( |26] |. Next, we 
describe a procedure to estimate the amplitude. 

In general, amplitude estimation is an important task at 
the receiver that is useful for interference cancellation and 
power control. The proposed blind interference estimation and 
tracking algorithm needs some form of amplitude estimation 
in order to accurately compute the interference power To 
estimate the amplitudes of the associated user signals, we 
describe the following procedure to estimate the absolute value 
of the output of the RAKE receiver defined in ( l30b as given 
by 

qk[l + l] = {l- P)qk[l] + \xk[A\- (41) 

The amplitude can be estimated by removing the square-root 
of the interference power from the above estimate according 
to 

Ak[i + !] = (!- /3)ifeM + (kfeMI - \vk[m (42) 

The above procedure is simple and effective to estimate the 
amplitude for use in the interference power estimation pro- 
cedure. Since one recursion depends on the other, a designer 
shall start the procedure with an interference power equal to 
zero (or equivalently \vk[i] \ — 0). 
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V. Analysis of The Algorithms 

In this section, we study of the properties of the opti- 
mization problems associated with the design of the SM- 
CCM-based algorithms and examine the convergence and 
tracking performances of the proposed algorithms. In order 
to ensure the convergence of the SM-CCM-RLS a persistence 
of excitation condition on the received data r[i] must hold 
and the transmitted symbols bk[i] and the noise n[i] have 
to be uncorrelated. To this end, the variable bound must 
adapt such that the above conditions are met, following the 
same procedure estabhshed for OBE algorithms in llTSl . Since 
the SM-CCM-RLS algorithm is expected to converge without 
misadjustment to the Wiener filter and certain aspects of the 
convergence analysis of the SM-CCM-SG algorithm are of 
greater interest, we will focus on the latter The CCM-based 
algorithms are inherently nonlinear and deal with time-varying 
environments, which leads to difficulties in the study of their 
performance. For this reason, we will resort to an effective 
approach termed energy conservation principle 081 that lends 
itself to such analysis. 

A. Analysis of the Optimization Problem 

Let us consider the optimization problem that needs to be 
solved for the design of the blind receiver and can be solved 
by the algorithms proposed in Section III: 

minimize JcM{wk[i]) = E[{w^ [i]r[i]\^ - if] 

subject to = u 

whenever - l)^ > 

(43) 

Let us rewrite the received vector as 

K 

r[^=Y. Akbk[^Pk\i] +J7fc [«] + n [»] , (44) 

where p^[i\ = Ckhk[i\, R = S + G + a"^!, S = 
E[Y.k=i M^^mi and G = i?Ef=i [z]], are 

independent and identically distributed random variables, and 
are statistically independent from n\i\. 

Consider user k as the desired user and let Wk[i] be the 
receive filter for this user Our strategy to analyze the above 
optimization problem and its properties is to transform the 
variables and rewrite the problem in a convenient form that 
provides more insight about the nature of the problem. Let us 
now define the signal of the desired and the signal of all the 
users after applying the receive filter 

tk^AkP^Wk[i], and t^AP''wk = [ti...tKf, (45) 

where P = [pi...pj^], A = diag[Ai . . . and b = 
[bi . . . 6a']^- The relation between the receive filter, the chan- 
nel and the signature sequence can be written as 

w^[i]Pk[i]^w^[i]Ckhk[i]^,y. (46) 



Therefore, we have for the desired user k the following 
equivalences 

tk[i] = Akpf!mwkm = Akh^ii] Cf t^feW = Afc/if 

hkli] 

(47) 

By considering the noise and the ISI negligible we can write 
the cost function in ( l43T l as 

JcmW = £;[kfeMI^-2|zfeH|2 + i] 

= E[{t"bb"tf] - 2E[{t"bb"t\ + 1 

i=i j=i i=i 

K K (48) 

K 

where t = f{wk[i]) is a linear function of the receive filter 
Wk[i], the terms multiplying the summations in the third line 
of (081) are obtained by evaluating the expected values in the 
second line of (gUl lUll, and D = tjt* = iy'^\Ak\'^\h'i! \i]hk[i\\. 
The strategy we employ to enforce the convexity of the 
optimization problem relies on the adjustment of the parameter 
1/ and a transformation of variables on the cost function in (l48T l, 
which will be detailed in what follows. 

Let us now transform the above cost function taking into 
account the constraint = ly and the fact that 

we are interested in demodulating user k and rejecting the 
remaining users. Therefore, we introduce another parameter 
vector i ~ [ti . . . tk-i ife+i . . . t^]"^ = Bwk[i] that excludes 
the user k and is responsible for the interference suppres- 
sion task of the remaining users (all but user k), where 
B = AP", P[i] = [piH . . . p^+iH . . .p^W] and 

A = diag[^i . . . Ak-i A^+i ■ ■ ■ Ak]- The transformed cost 
function is given by 

JcM{t) = KD + i"tf - AiD^ + t"tf ^^^^ 
-'i{D + t"t) + l. 

At this point we need to take into account the constraint 
JcM{wk[i]) < 7fe[*]- Since i = Bwk[i] is a linear mapping 
we have an equivalent constraint JcAiii) < 7^[i], where 7^[i] 
is the bound modified by the linear mapping. It can be verified 
that the constraint set generated by JcAiii) < Iki^] is not 
convex and leads to two disjoint parallel hyperstrips in the 
parameter space. 

Given that the constraint is not convex, the optimization 
problem is clearly non convex. However, in this context non- 
convexity only poses a problem if local minima are present 
and prevent an algorithm from reaching the global minimum. 
It turns out that a designer can adjust the parameter u in order 
to enforce the convexity for each hyperstrip. Computing the 
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Hessian ||43] H 



dt 



dt 



H = 1&[{D- 1/4)/- 



t"tl 



-diag(|tiHp...|tfc_i|2 \tk+i\ 



positive or positive semi— definite terms 



The condition D = v'^\Ak\'^\ 



h"[i\hk[\ 



(50) 

P > 1/4 ensures that 



where R^^ = £:[efc[i]t>fe[i]r^[i]] and R^rWo^t = 0. 

From the above, it can be concluded that Wk[i] converges 
to Wk.opt and ( fSSl ) is stable if and only if ni^o(-^ ~ 
]i?[^fc[z]]ili,r) — 0, which is a necessary and sufficient 

'condition for lim^^oo £^[eiu J = and Wk.opt- 

For stability, a sufficient condition for (ISSl l to hold implies 
that 

0<£;[//(oo)] <min^, (56) 



. I A- 1' 

where A^.'' is the fcth eigenvalue of R^r that is not necessarily 
real since R^, is not symmetric. 



there is no local minimum in each of the hyperstrips resulting 
from the optimization problem and convexity can be enforced 
in each of the hyperstrips. Since t = Bwk[i\ is a linear 
function of Wk\i\ then JcM{i) preserves the maxima and 
minima properties of JcMi^kii])- For sufficiently high signal- 
to-noise ratio (SNR) values, the extrema of the cost function 
can be considered a small perturbation of the noise-free case 

m- 

B. Stability Analysis 

In this part, we discuss the stability analysis of the SM- 
CCM-SG algorithm described in subsection III.A. In particu- 
lar, we consider the range of step-size values for convergence. 
Let us now rewrite the update equation of the algorithm as 

Wk[i + 1] = Uk[i\{wk[i] - iik[i\ek[i]zl[i\r[i]) + vPk[i]{Pk[i\PkiMASE receiver Consider the MSB at'time i: 



C. Steady-State Analysis 

In this part of the analysis we are interested in devising 
a formula to predict the excess MSB, which depends on the 
MAI, the ISI, the noise at the receiver and the parameters of 
the SM-CCM-SG algorithm. The excess MSB is related to the 
error in the filter coefficients e^j^, [i\ via the a priori estimation 
error, which is defined as 



(57) 



where [i] — Wopt — Wk[i], and lOopt is the optimum linear 



= Wk [i] - fik \i]ek [i]r" \i]wk [«](/- 



Pk\i]Pk\i] 



ir z 



MSE[i] 



= Wk[i\ - ^ikmekmv[i\r [i\Wk[i\ 
= {I - fik[i\ek[i\v[i\r"[{\)wk[i\, 



E[\bk[^^^k 

e^in + E[\ea[i]\'']+p^\i]E[ 



(51) 



i]]+E[e^^[i]]Pk[i] 

A^]~E[eZ,[i]r[i]r''[i]wk,opt], 
(58) 



where Vk[i] = (/ — ^h^I^jJj '''[*]) ^^'^ the expression in 
the second line of (ISTI i is obtained by substituting IIk[i] = 
/ — Pk[i]{Pk[i]Pk[i])^^P^[i\ into the first line and further 
manipulating the terms. Let us now define the error vector 

e-ujj* + 1] = iffe,opt - Wk[i + 1] 



where 



^opt' 



and £'[e„Ji]] 



(wopt - Wk[i])"r[i] and e,nin = 
When i oo, we have if — !■ Wk^opt 
and the steady-state MSB 



lim MSE[ 



lim E[\ea[i]f] 



fJ.k[i\ek[i\Vk[i\r"[i])ei^J{\ + fik\i]ek[i\vk[i\r" [i\Wk.opt 

(52) 



The steady-state excess MSB is then defined [I] as 



lim E[\eam' 



(59) 



(60) 



In order to proceed with the analysis, we need to resort to an 
assumption. 

Assumption 1: Let us suppose that for the algorithm in (ISTj 
when i — > oo 



E[fXkm[q]^E[^ikm[4[^ 



(53) 



This assumption holds if fik [i] is a constant, and we claim that 
it is approximately true if fik [i] varies slowly around its mean 
value. By writing 

= E[^ikm[elll]] + E[U,k[i] - i?[/ifeM])e|M], 

(54) 

we can notice that the second term on the right-hand side will 
be small compared with the first one provided that ^k [i] varies 
slowly around its mean value. 

By taking expectations on both sides of (l52b and using the 
previous assumption, we have 



E[e, 



1]] = {I - E[^ik[i]]Rrv)E[e, 



(55) 



Using the energy conservation principle 11381 . the proposed 
SM-CCM-SG algorithm can be written in the form 



Wk[i+1] = Wk[i] 



Pk[i]Pk[i] 



(61) 

where Fg^, [i] is a generic scalar function determined by the 
adaptive algorithm. Subtracting the above recursion for Wk[i + 
1] from Wk. opt, we obtain 



\i + l] = 



- Hk[i\uk[i\Fe^[i\ 



Using the a priori estimation error ea[i] = 
ing the previous equation, we obtain 



ir I 



ea[i]^ep[i]+fik[i]uj^[i]r[i]F*[ 



ep[i\+ ^lk\i]r"\i]{I 



(62) 
Rewrit- 

(63) 
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Since r^\^ I-EMpj^ 



r[i] = uj! [i]uk[i], we have 



fik\i]u"[i\uk\i]F* 



(64) 



Using the fact that F*^[i] 
e,„Ji] as 



we can rewrite 



into ( |7Q] i. we obtain 

+ 3E[^l[i]]E[uj![i]uk[i]]Em^']E[Km'] 
+ 3E[f^l[m^k[i]nk[i]]E[Am[\eam'] 

+ E[^'[l]]E[u^\^]u,mm^^'m + m^m 
E[^^imK[^]ukm+E[f^im[un^Mm[Km 



e^, [i + l] = e^fc [i] - 



u"\i]uk\i] 



(65) 



2i;[MfcW](S[MAl2]i;[|e„Wp]+i;[i,2[,]]s[|e„[i] 
+i?[r,2W]i?[|e,Hp]+i?[|e4z]|4], 



Rearranging, we obtain 



D 



(66) 



where A 



Scr'^MAI 

— ^6 



Zat and B = a! 



'MAI 



"^^MAI^ti 



+ ^MAI + '''mAI 



2ct 



MAI"ti 



^MAI 



4a, 



By defining = l/(it^[i]itA;[i]), squaring the previous 
equation and taking expectations, we obtain the relation 



2al 



we can 

JlmA 



(71) 
. -6^2 + 

^O^MAI + 2- 

that 



assume 

,2n 



i?[||e„,[z+l]||^]+i?[/2fe[z]|epH|^] 
In the steady state, we can write 

,M\?]=E[\\e^A^+m\ 



= E\ 



N+1 



(67) 



Upon convergence(i — > oo), 

Cr^™, and = (£;[w2])m ^ ^2m^ ^j^gj-g ^^^^^^ ^^^^ 

and^i, are the variances of the Gaussian distribution. In this 
'gikuatibn, the high power terms £'[|eci[i]|"'] and £'[|eci[i]|^] may 
be neglected as these values are typically very small compared 
with the remaining terms. The excess MSE is obtained as 
follows 



E 



and E[liMeM^] = E[fjLk[i]\eM?] 

(68) S, = E[\ea[if\ 



i?[M?.(oo)]i?[<W«fc[j]]S 



The energy preserving equation can be used to compute 
the excess MSE. It can be obtained by cancelling the terms 



2E[iik{'x)]{G 



1 

MAI 



al^ol)-E{i,l[^)\E{u^{i\uu 
(72) 



i?[||e^,[z + l]||2]andii;[||e^J 
(|64] | into (l63, which yields 



in Wi\ and by substituting 



E{M^\ 



= E\lLu\^ 



e:a^ 



(69) 



where Ff,^. [z] 



{\~Wm^)zl{i\ and zM 



MAI + ri^.[i] + v[i] - ea[i]- Substituting Fe^^l 
and manipulating the terms, we obtain 



Akbk\i] + 
and ( l64b into 



+E[pikmw[m[ea[^4mi- 



\ZkW 



E[^,l[l]]E[u^[^ 



c 

Uk[i\\zk[i] 



(70) 



\zk m 



At this point we need to resort to another assumption to 
continue with our analysis. 

Assumption 2: In the steady state, tt^[i]Mfe[i] and li^e^WP 
are uncorrected. The quantities MAI, rj[i], ea[j]} 

are zero mean random variables, and are mutually independent. 
We use the fact that E[b1"'] = 1 for any positive integer m and 
that the residual MAI and ISI are Gaussian random variables. 

Using the previous assumption and substituting Zk[i] ~ 



D. Tracking Analysis 

In this subsection, we assess the proposed SM-CCM-SG 
algorithm in a non-stationary environment, in which the algo- 
rithm has to track the minimum point of the error-performance 
surface. Specifically, we derive an expression for the excess 
MSE of a blind adaptive linear receiver when the channel 
varies in time. Differently from the works in ll38l . Il40l . where 
expressions were derived for the constant modulus (CM) 
algorithm and the CCM algorithm, respectively, we consider 
a set-membership approach. In the time-varying scenarios of 
interest, the optimum receive filter coefficients are assumed to 
vary according to the model tUopt + 1] = '"'opt + where 
q[i] denotes a random perturbation. This is consistent with 
tracking analyses of adaptive filtering algorithms and requires 
an assumption. 

Assumption 3: The sequence q[i] is a stationary sequence 
of independent zero-mean vectors and positive definite auto- 
correlation matrix Q ~ E[q[i]q^ [i]], which is mutually inde- 
pendent of the sequences {^^[z]}, {MAI[i]}, {77[i]}. 

Let us consider the weight-error vector e^^ [i] = Wk.opt — 
Wk[i], which satisfies 

£-Wk + M - Mfe W Wfe \i]Fe^ [i] +q[i\. (73) 

With ea[i] = e^fc W^W ™d SpW = ^Si-I* + we have 



(74) 
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Using ( TTST l and d74] i. we obtain 

+ 1] + fik[i]uk[i]el[i] = e™J«] + /2fc[i]Mfc[i]e*[z]. 

(75) 

Squaring (|75] l and taking the expected value on both sides, we 
obtain 



E[fl,[i]\ea[i]\^]^E[\\e^,[^+qm 
+ E[j2k[i]\ep[i]\% 



(76) 
where 

i;[||e„,[z] + q[z]||2] = i?[||e„,[z]||2] 

+ £;[q«We„,W]+£;[q^[%[z]]. 

(77) 

Using Assumption 3, we have £^[e:^^ = 
E[q"[i]e^,[i] = 0. When ^ ^ oo E[\\e^,[t + = 
-^[Ikiufc '^he energy preserving equation describing the 

tracking performance is given by 

i?[/2fcW|e.W|2] = Tr(g)+i?[/2,.W|eaW-^F;jz]p (78) 

Expanding the equation above, it can be simplified to 

C = Tr(Q) + D, (79) 

where C and D are defined in (fTTT i. The excess MSE is then 
obtained as 



E[^il{^)]E[\\uk[iW]B + TiiQ) 



2i?[/x,(cx))KAj + a2 + a2) - i?[^2(cx3)]i?[|w,[z]P]A 



(80) 



£. Computation of Moments 

In order to compute the expressions obtained in (l72t and 
dSOl l, we need to obtain the first and second-order moments of 

[oo] ■ To this end, we resort to the expression for the variable 
step size in ( |20] | and the methodology employed in |l3l, which 
after some algebraic manipulations yields 

(1-Pup) E[Ak[^] 



S[/ifc[oo]] =S[7fc[i]]Pup 



E[f,l[^]] ^ E[jk[i]]Pup 



E[^M E[\u,m'] 

(1-Pup) E[Al[q] 



where the probability of update Pup is given by 

Pup = Pr[\ek\i]\ > P[7fcW]] = Pr[\ek\i]\'' > E[\lk[i'. 



(81) 
(82) 



= 2Q 



(83) 



where Pr[-] denotes the probability Q{x) is the complemen- 
tary Gaussian cumulative distribution function P4l given by 



Q{x) 



12-K 



(84) 



VI. Simulation Results 

In this section we assess the performance of the pro- 
posed and analyzed adaptive algorithms in terms of mean- 
square error (MSE) and bit error rate (BER). In particular, 
we consider the proposed SM-CCM-SG and SM-CCM-RLS 
adaptive algorithms and the existing CMV-SG and CMV- 
RLS reported in ||26l, the CCM-SG 1321 and the CCM- 
RLS Ea algorithms, the SM-CMV-SG and SM-CMV-RLS 
reported in JT] with and without the proposed blind PDB and 
PIDB time-varying bounds described in Section IV. The DS- 
CDMA network employs Gold sequences of length N = 31, 
the users are randomly distributed and communicate over 
multipath channels. The channels experienced by different 
users are independent and different since we focus on an 
uplink scenario. The DS-CDMA system under consideration 
employs quadrature phase-shift keying (QPSK) modulation. 
The channel coefficients are given by hkj[i] = Pk.iO!k.i[i], 
where ak,i[i] {I = 0, 1, . . . , — 1, k ~ 1,2, . . . , K) are 
obtained with Clarke's model B4l and pk.i represent the 
powers of each channel path. We show the results in terms 
of the normalized Doppler frequency f^T (cycles/symbol), 
where fd is the maximum Doppler shift and T is the symbol 
interval. We use three-path channels with relative powers given 
by 0, —3 and —6 dB, where in each run the spacing between 
paths is obtained from a discrete uniform random variable 
between 1 and 2 chips. The proposed blind channel estimator 
described in Section IV and that of II41I model the channel 
as a finite impulse response (FIR) filter and we employ a 
filter with 6 taps as an upper bound for the experiments. The 
phase ambiguity derived from the blind channel estimation 
method in BTI is addressed in our simulations by using the 
phase of hfc(O) as a reference to remove the ambiguity and 
for time-varying channels we assume ideal phase tracking. 
Alternatively, differential modulation can be used to account 
for the phase rotations. The tuning parameter a, the forgetting 
factor /3 and the weighting parameter r required by the time- 
varying bounds described in Section IV have been obtained 
by experimentation and chosen such that the performance 
of the algorithms is optimized. The update rate (UR) has 
been computed by counting, for each simulation trial t, the 
number of updates (Nu^t) performed and then dividing it by 
the number of received symbols ( A^s,*)- Then, the URis given 
hy UR= ^ J2l where T 1000 is the total number 
of trials. 



A. MSE Analytical Results 

The aim of this part is to verify the validity of the analytical 
results obtained in Section V. Specifically, we shall evaluate 
the analytical MSE obtained with the analytical formulas in 
(|72l) and and compare them to the results obtained by 
simulations. We consider first a scenario with fixed channels 
ifd = 0) and assess the MSE curves using ( |72] | against the 
number of received symbols and also versus the signal-to- 
noise ratio (SNR) defined by Ei^/Nq, as shown in Fig. |2] The 
results show that the analytical curves match very well those 
obtained via simulations, showing the validity of our analysis 
and assumptions. 
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b) N=31,K=10 users 
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a) N=31, EJN„ = 12dB, K=10, f T=5x10 
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Fig. 2. MSB performance against a) tlie number of received symbols at 
Ef,/No = 15dB and b) Ef,/No (dB). 



Fig. 4. Performance of the a) interference power estimation and tracking 
and b) amplitude estimation algorithms at Ei,/Nq = l2dB. 



The tracking analysis of the proposed SM-CCM-SG algo- 
rithm in a time-varying multipath fading channel is discussed 
next. We consider the same scenario as before except for the 
fact that the channel is now time-varying and the analytical 
results of the tracking analysis in (ISOl l are employed. We 
consider a Jakes' model with a typical normalized fading rate 
fdT. We compute Tr(Q) with the aid of Jo{2iifdT) B4l . 
which is the autocorrelation function of the Jakes' model and 
where Jq is the zero-order Bessel function of the first kind. A 
comparison of the curves obtained by simulations and by the 
analytical formulas is shown in Fig. [3] Similarly to the case 
of fixed channels, the analytical and simulated curves agree 
and show the validity of the proposed formulas. Nevertheless, 
the curves indicate a higher misadjustment as compared to the 
curves in Fig. [2] due to the time-varying process and extra 
effort of the adaptive filters to track the channel variations. 



a) N=31, K=10 users, E /N =15dB, fj=10 



b) N=31, K=10 users 



-^i>— Simulated(a=5,T=0.25, p=0.£ 
— Analytical 
- Optimal 
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Fig. 3. MSE perfomiance against a) the number of received symbols at 
El, /No = IbdB and b) Ei,/No (dB). 



B. Interference Power and Channel Estimation 

At this point, we wish to evaluate the effectiveness of 
the proposed algorithms for estimating and tracking the in- 
terference power and the amplitude. To this end, we have 
carried out an experiment, depicted in Fig. |4] where the 
proposed algorithm estimates of the MAI and ISI powers 
and the amplitude of the desired user are compared to the 
actual interference power and amplitude. The results show that 
the proposed algorithms are very effective for estimating and 
tracking the interference power in dynamic environments, as 
depicted in Fig. |4] 

We assess the proposed channel estimation (CE) algorithm, 
called SM-CCM-CE, with the time-varying bounds PDB and 
PIDB, and also with a fixed bound, and compare them to 
the CMV-based method (CMV-CE) reported in lEI and the 
subspace algorithm of ll42l in terms of the MSE between 
the actual and the estimated channels using the following 
dynamic scenario. The system has initially 10 users, the 
power distribution among the interferers follows a log-normal 
distribution with associated standard deviation of 3 dB. After 
1500 symbols, 6 additional users enter the system and the 
power distribution among interferers is loosen to 6 dB. The 
results, shown in Fig |5] reveal that the proposed SM-CCM- 
CE algorithm outperforms the CMV-CE technique reported in 
ll4n because SM-CCM-CE uses a variable forgetting factor. 
The proposed SM-CCM is also computationally simpler than 
CMV-CE due to the sparse updates, whereas it is substantially 
less complex than the subspace method of |42|. This is because 
the SM-CCM-CE needs to compute the principal eigenvector 
of the Lp X Lp matrix CkRk [i]Ck, while the subspace 
technique in ll42l requires the principal eigenvector of the 
]\I X M matrix R [i]. Since Lp ^ M in most practical 
scenarios SM-CCM-CE is typically considerably simpler than 
the subspace method. 

C. BER Performance 

The SM-CCM algorithms are assessed in a non-stationary 
environment where users enter and exit the system, as depicted 
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Fig. 5. MSB performance of channel estimation versus number of received 
symbols in a dynamic scenaiio where receivers operate at SNR = E^/Nq = 
15 dB for the desired user. 
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Fig. 6. BER performance of the algorithms versus number of symbols for 
a non-stationary scenaiio. The parameters are 7 = 1.3 for the SM-CMV, 
7 = 0.65 for the SM-CCM. 



in Figs. |6] and |7] The system starts with 2 interferers with 7 
dB above the desired user's power level and 5 interferers with 
the same power level of the desired one, which corresponds to 
the signal-to-noise ratio Eh/No = 15 dB. At 1000 symbols, 
2 interferers with 10 dB above the desired signal power level 
and 1 with the same power level enter the system, whereas 
1 interferer with 7 dB above the desired signal power level 
leaves it. At 2000 symbols, 1 interferer with 10 dB above, 
and 5 interferers with the same power level of the desired 
signal exit the system, while 1 interferer with 15 dB above 
the desired user enters the system. The results for 100 runs 
show that the proposed SM-CCM-RLS algorithm achieves 
the best performance, followed by the proposed SM-CMV- 
RLS recursion, the CCM-SG and the CMV-SG methods. In 
summary the SM-CCM algorithms outperform the SM-CMV 
techniques in all scenarios and the SM-CCM-RLS algorithm 
is the best among the analyzed algorithms. 

In a near-far scenario the eigenvalue spread of the covari- 
ance matrix of the received vector r[i] is large, affecting the 
convergence performance of the SG algorithms with fixed step 
size and making it very difficult to compute a pre-determined 
value IT] for the step size. In this case, the SM-CCM algo- 
rithms are able to deal with near-far situations since they adopt 
variable step size or variable forgetting factor mechanisms, 
ensuring good tracking performance in dynamic scenarios 
and an improved performance over the existing CCM-type 
algorithms. In addition, due to their data-selective update 
feature the SM algorithms can save significant computational 
resources as they only require parameter updates for about 
20% of the time for the SG-type recursions and around 12% 
for the RLS-based techniques. 

The same scenario illustrated in Fig. |6] is considered for 
the SM-CCM algorithms with time-varying error bounds, as 
shown in Figs. [8] and |9] The results indicate that the proposed 
time-varying bounds are capable of improving the performance 
of SM algorithms, while further reducing the number of 
updates. The RLS-type algorithms outperform the SG-based 
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Fig. 7. BER performance of the algorithms versus number of symbols for 
a non-stationary scenaiio. The parameters are 7 = 1.3 for the SM-CMV, 
7 = 0.65 for the SM-CCM. 



techniques as expected and have lower UR. Moreover, the 
algorithms with the PIDB approach achieve the best perfor- 
mance, followed by the algorithms with the PDB method and 
the SM recursions with fixed bounds. With respect to the UR, 
the PIDB technique results in the smallest number of updates, 
followed by the PDB approach and the fixed bound. 

The BER performance versus Ei,/No and number of users 
is shown in Fig. [TOl We consider data packets of P = 1500 
symbols and compare the proposed SM-CCM-RLS algorithm 
with a fixed bound, with the PIDB mechanism, the training- 
based BEACON algorithm with the PIDB mechanism ||2l] and 
the linear MMSE receiver that assumes perfect knowledge of 
the channels and the noise variance. We have measured the 
BER after 200 independent transmissions and the BEACON 
algorithm 161. 121] is trained with 200 symbols and is then 
switched to decision-directed mode. The curves illustrate that 
the proposed bUnd SM-CCM-RLS algorithm can approach the 
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Fig. 10. BER performance versus (a) E^,/Nq with K = 10 users and (b) 
Fig. 8. BER performance of the SG algorithms versus number of symbols number of users (K) at Et,/No = 15 dB . 



for a non-stationaiy scenario with time-varying bounds. The parameters are 
a = 8, T = 0.35 and /3 = 0.95 for the time-varying bounds. 
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Fig. 9. BER performance of the RLS algorithms versus number of symbols 
for a non-stationary scenario with time-varying bounds. The parameters are 
a = 8, T = 0.35 and /3 = 0.95 for the time-varying bounds. 



performance of the supervised BEACON algorithm and that 
of the linear MMSE receiver 

VII. Conclusions 

We have proposed SM-CCM adaptive algorithms for inter- 
ference suppression in DS-CDMA systems. We have analyzed 
the optimization problem that gives rise to the SM-CCM 
algorithms and devised analytical expressions to predict their 
MSE performance with a good accuracy. The proposed SM- 
CCM algorithms are equipped with variable step sizes and 
forgetting factors and are only updated for a small fraction of 
time without incurring any performance degradation. A blind 
framework for SM estimation that takes into account parameter 
estimation dependency and MAI and ISI for multiuser commu- 
nications has also been introduced to provide a time-varying 
bound that is robust against channel and SNR variations. Sim- 
ulations have shown that the proposed SM-CCM algorithms 
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